Configuration mixing of angular-momentum and particle-number projected 
triaxial HFB states using the Skyrme energy density functional 
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We present a method based on mean-field states generated by triaxial quadrupole constraints 
that are projected on particle number and angular momentum and mixed by the generator coordi- 
nate method on the quadrupole moment. This method is equivalent to a seven-dimensional GCM 
calculation, mixing all five degrees of freedom of the quadrupole operator and the gauge angles for 
protons and neutrons. A first application to '^''Mg permits a detailed analysis of effects of triaxial 
deformations and of K mixing. 



I. INTRODUCTION 

Methods based on the self-consistent mean-field ap- 
proach [l[ are up to now the only microscopic tools that 
can be applied to all nuclei including the heaviest ones. 
However, mean-field methods have several limitations, 
and the method that we present here is part of an at- 
tempt to eliminate two of the most penalizing presented 
in a series of papers . The first limitation is due to the 
determination of a wave function in an intrinsic frame 
of reference. Although the symmetry-breaking mean- 
field approach is a very efficient and transparent way to 
grasp the efi^ect of correlations associated with collective 
modes in the limit of strong correlations [1, d, 0, Q , the 
absence of good quantum numbers and the correspond- 
ing selection rules does not allow direct determination 
of the electromagnetic transition probabilities. Instead, 
approximations have to derived based on the collective 
model and these cannot cover all possible cases as they 
are well justified only in the limit of large deformation. 
The second limitation concerns nuclei for which a mean- 
field description through a single configuration breaks 
down because several configurations with different shell 
structure are close in energy without being separated by 
a large potential barrier. Methods to overcome these 
two problems have been proposed in the past, but it is 
only in the last ten years that sufficient computational 
resources have become available to construct and apply 
methods based on realistic effective energy density func- 
tionals aiming at eliminating both of these limitations. 

Angular momentum projection 0| is the key tool 
to transform the mean-field wave function from the in- 
trinsic to the laboratory frame of reference. There is no 
ambiguity in determining directly electromagnetic tran- 
sition probabilities when working in the laboratory sys- 
tem. However, without the simplifying assumption of ax- 
ial states, the restoration of rotational symmetry requires 



that rotations about three Euler angles be considered. So 
far, for mean-field states with triaxial quadrupolar defor- 
mations, this has been mostly done in the context of phe- 
nomenological models by using small shell model spaces 
and often schematic interactions 

the context of energy density functional methods, pro- 
jection on angular momentum of mean-field states with 
triaxial quadrupolar deformations has been performed 
at the Hartree-Fock level by using a simplified Skyrme 
interaction in Ref. [T^ and very recently with a full 
Skyrme energy functional in Ref. p^ . Cranked wave 
functions have been projected in both cases to approx- 
imate a variation after projection on the angular mo- 
mentum [l^ . [tR [Tsj . but pairing correlations were not 
included. 

To solve the problem of energy surfaces that are soft 
with respect to a collective degree of freedom, one has to 
introduce fluctuations in this collective degree of freedom 
into the ground-state wave function using the generator 
coordinate method (GCM) H, i, 0, [H. A gam, many 
studies used small shell model spaces and schematic in- 
teractions, but there is also a quite a large body of work 
starting from mean-field methods based on effective inter- 
actions in the full model space of occupied single-particle 
states and with inclusion of pairing correlations. There 
have been several applications [2^ l2l|, that deal with 
the intrinsic quadrupole mode including triaxial defor- 
mations, which is in most cases the dominant low-lying 
collective excitation mode in nuclei. There also have been 
studies of octupole mod es l23l[2^ [25| and their coupling 
to the quadrupole mode [2^, [23| , as well as investigations 
of fluctuations in pairing degrees of freedom [H, [2^, [s^] ■ 
All the studies mentioned so far have in common the fact 
that simultaneous symmetry restoration, if performed at 
all, is limited to numerically inexpensive modes, such as 
particle number or parity, whereas angular momentum is 
not restored. 

A simultaneous treatment of angular momentum pro- 
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jection and of fluctuations with respect to triaxial 
quadrupole deformation requires that five degrees of free- 
dom be considered. The most transparent representation 
consists of three Euler angles defining the relative orien- 
tation of the intrinsic major axis frame in the laboratory 
frame and two independent degrees of freedom character- 
izing the intrinsic deformation (e.g., through the coordi- 
nates (3 and 7 introduced by A. Bohr). In the context 
of the GCM, however, this has so far only been done in 
highly schematic models, by using either a single-j shell 
3ll | or a very small shell model space and a schematic 
interaction [s^. [ssj. 

A first step toward the simultaneous treatment of 
these five degrees of freedom has been carried out re- 
cently with the GCM mixing of angular-momentum pro- 
jected quadrupole-deformed axial mean-field states. This 
scheme is nowadays routinely applied by several groups, 
using Skyrme energy density functionals_ S^, Issl [sg'I , the 
density-dependent Gogny forc e .37i , .3S , .39j , or a relativis- 
tic point-coupling model [13, EfT In all of those cases, 
intrinsic triaxiality is neglected and two of the three Euler 
angles can be treated analytically. 

All the studies mentioned so far that consider angu- 
lar momentum projection have in common the fact that 
this operation is performed after variation. This has 
several drawbacks Q, in particular when working with 
time-reversal invariant mean-field states. The alterna- 
tive, variation after projection, is computationally very 
demanding. The only standard method where variation 
after angular-momentum projection is considered (to- 
gether with restoation of A^, Z, and parity) is the VAM- 
PIR/MONSTER approach [H, El] which is confined to 
a very small shell model space and a shell model effective 
Hamiltonian. In this framework, the intrinsic dynamics 
is not described in terms of deformation degrees of free- 
dom, but rather in terms of quasiparticle excitations. 

An approximate scheme to describe five-dimensional 
quadrupole dynamics in a fully microscopic framework 
has been set up long ago. The idea is to construct the in- 
put of the Bohr-Hamiltonian microscopically j^, [l^, [3] . 
This is an approximation to the full five-dimensional 
GCM as it replaces the nonlocal kernels entering the Hill- 
Wheeler-GrifRn equation for nonorthogonal weight func- 
tions by local potentials and mass parameters in a collec- 
tive Schrodinger equation for orthogonal collective wave 
functions, but it also allows one to incorporate improved 
moments of inertia at a moderate cost. This approach 
has routinely been used for a long time in the framework 
of the density-dependent Gogny force [i^, \^ , and it has 
also been recently set up for Skyrme interactions [47j . 

In this paper we present a first step toward a micro- 
scopic treatment of triaxial quadrupole dynamics and of 
angular momentum projection in the context of nuclear 
energy density functional methods using the full space 
of occupied single-particle states. Th emethod employed 
generalizes that described in Refs. [13, [H, [s^ by allow- 
ing for the breaking of axial symmetry of the mean-field 
wave functions. 



II. THE METHOD 

For a given nucleus, the calculations are performed 
in three steps. First, a nonorthogonal basis of self- 
consistent mean-field states is generated with constrained 
Hartree-Fock-Bogoliubov (HFB) calculations. Second, 
angular-momentum and particle-number projected ma- 
trix elements between all mean-field states are calculated. 
Third, these matrix elements are used in a configuration- 
mixing calculation to determine the correlated ground 
state, the spectrum of excited states, and the transition 
moments 

The only inputs to the calculation are the proton and 
neutron number and the parameters of a universal energy 
functional. The latter are taken from the literature and 
obtained from a global fit to nuclear properties aiming at 
the description of the entire chart of nuclei, without local 
fine tuning. One consistency requirement of our method 
is that the same effective interaction is used to generate 
the mean-field states and to calculate the configuration 
mixing. We chose a Skyrme energy functional supple- 
mented by the Coulomb interaction in the mean-field 
channel, together with a zero-range, density-dependent 
functional in the pairing channel. 

We will give here, and in the Appendix, only a sketch 
of the ingredients of the method. We postpone a detailed 
description to a future paper on the generalization of the 
code to time-reversal invariance breaking triaxial states. 



A. The mean-field basis 

The HFB equations are solved self-consistently with 
the two-basis method of Ref. . The quadrupole mo- 
ment is constrained through two coordinates qi and (72, 
which are related to the usual mass quadrupole defor- 
mations (3 and 7 of the Bohr Hamiltonian through the 
relations 14911: 



91 
92 



Qo cos(7) - Qo sin(7) 



with 



V3 



P2 



h sin(7) , 



5 AnQo 
16n 3R^A 



(1) 
(2) 

(3) 



where R = \.2A^^^ fm. A mesh of positive values 
for q\ and qi covers the entire first sextant of the /3-7 
plane. The mean-field equations are solved on a three- 
dimensional mesh [4^, with the total nuclear density 
imposed to be symmetric with respect to three planes. 
Thus, all odd-Z moments of the density identically vanish 
whereas the even moments, which are not constrained, 
are fully taken into account and take the values that min- 
imize the energy. All the wave functions that we consider 
are, therefore, restricted to have positive parity. 
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For the results described in the following, the mean- 
field wave functions are required to be time-reversal in- 
variant, such that the collective coordinates can be lim- 
ited to one sextant of the /3-7 plane. The five other sex- 
tants correspond to one or several permutations of the 
principle axes of inertia. After projection and configu- 
ration mixing, all of them are equivalent, although in- 
termediate (nonobservable) quantities might differ. This 
equivalence will be used as a test of the numerical imple- 
mentation of our method of Sec. III. 



B. Projection 

Restoring symmetries allows to extract states with 
good quantum numbers from the mean-field states and 
provides a clean framework to use selection rules for 
electromagnetic transitions. Eigenstates of the particle- 
number operator are obtained in the same way as pre- 
viously [3^, [3^. From a technical point of view, it is 
crucial to restore this symmetry when mixing wave func- 
tions that are not eigenstates of the particle-number op- 
erators. Since those wave functions have only the right 
mean particle number, their transition matrix elements 
will contain nondiagonal contributions in the number of 
particles. One can estimate easily the error that this sym- 
metry breaking brings to energies. The binding energy 
per nucleon of most nuclei is between 7 and 8 MeV/u 
and a deviation of the particle number of a mixed state 
as small as 0.1 nucleons would affect its energy by several 
hundreds of keV. Particle-number projection removes 
this ambiguity (2^ . A technically simpler approximate 
method would be to introduce a constraint on the aver- 
age value of the particle number Nq of the mixed states 
in the Hill- Wheeler- Griffin equation [2^ [3^, but since 
the same value for this constraint has to be used for all 
the eigenstates to maintain their orthogonality, the error 
on the mean number of particles will vary from one state 
to another. We always project on the proton and neu- 
tron numbers imposed on the mean-field wave functions 
by Lagrange multipliers in the mean-field equations. We, 
therefore, drop any reference to particle numbers for the 
sake of simple notation. 

Angular momentum projection is significantly more 
complicated when triaxial mean-field states instead of ax- 
ial ones are considered. Eigenstates of the total angular 
momentum operator in the laboratory frame, J^, and its 
projection on the z axis, Jz, with eigenvalues h^J{J + 1) 
and hM, respectively, are obtained by applying the op- 
erator 



and ■Dj(^j^(a, /9, 7) a Wigner function ^ [s^. Pjlfj^ picks 
the component with angular momentum projection K 
along the intrinsic z axis. The projected state is then ob- 
tained by summing over all K components of the mean- 
field state \q), 

+J 

\JMnq) = fi(K) Pi,i,\q) 

K=-J 
+ ./ 

= ^ fi{K) \JMKq) (5) 

with weights f^{K) determined by minimizing the en- 
ergy 0, [il (see the following) . The number of values that 
the index k can take is restricted by the symmetries of 
the mean-field states [11] (signature with respect to x 
and time-reversal invariance) to J -I- 1 for even J values 
and J — 1 for odd J values. 

Note that P^m '^o^ ^ projector in the strict math- 
ematical sense, but it has the properties 
5 J J' 5 MM' P'kk' and {Pkm)"^ 



pj 



51|. 



3 J pj' 

KM^M'K' 



C. Mixing of states with different deformations 

The fluctuations of the intrinsic deformation and the 
resulting spreading of the nuclear states in the (3-^ 
plane can be described by a superposition of angular- 
momentum projected states using the two intrinsic 
quadrupole degrees of freedom q as generator coordi- 
nates. Taking into account that one has also to mix all 
K values for each deformation, one obtains a resulting 
wave function given by 



\JMv) ^J2T. 9) \JMKq) 



(6) 



K 



This corresponds to the discretized version of the GCM. 

{K, q) is a weight function of the K components of 
the angular-momentum projected states of intrinsic de- 
formation (7, determined from 



{JMv\H\JMv) 



SFJ*{K,q) {JMv\JMv) 



= 0, 



(7) 



which leads to the Hill- Wheeler-Griffin (HWG) equation 

[Ulil 

q' K' 

-EiXj{K, q; K\ q')] {K\ q') = (8) 



Pmk = da dp sin(/3) Vi;^R, (4) 



Alternatively, the integration intervals can be chosen as a g 
[0, 27r], /3 S [0,7r], and 7 S [0, 47r]. In any case, for systems with 
. . . integer J values as discussed here, the integration over [2n, 4it] 

where R = e~'^°'''' e"*'^'^" e"'^"^^ is the rotation operator of either a or 7 gives just a factor of 1. 
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with the energy and norm kernels 

H\K,q-K',q') = {JMKq\H\JMK'q') (9) 
T\K,q;K',q') = {JMKq\JMK'q') . (10) 

For the sake of simple notation, we have introduced the 
method using a Hamilton operator H. We shall comment 
on the procedure used to calculate the Hamiltonian ker- 
nels from an energy density functional, and discuss some 
problems this may cause in Sec. IIIFI 

The dimension of the variational space is considerably 
increased by introducing triaxial deformation. One can- 
not limit the angular momentum projection to K ^ Q, 
and several K components for each deformation have to 
be mixed. This has to be done with some care because of 
the high redundancy of the GCM bases, in the K space 
for a given deformation on the one hand and also for 
the whole set of deformations on the other hand. For an 
efficient elimination of redundant states we have chosen 
a four-step procedure. We start from a basis \JMKq). 
First, we restrict the subspace for each deformation q and 
angular momentum J by diagonalizing the norm kernel 
I'^{K, q; K', g), rewritten in its significant subspace p^ . 
and neglecting eigenstates with negligible eigenvalues, 
typically smaller than 10"'^. Second, we solve the HWG 
equation in K space within this new basis. The solu- 
tions \JMKq) are labeled by an index k, Eq. (O. The 
number of values that k can take for a given J depends 
on the deformation. It is just one for axial states and 
even J values and is usually larger for triaxial deforma- 
tion. In a third step, we transform all matrix elements 
to the new basis \ JMKq) and diagonalize the norm ma- 
trix T"^ {k, q; k', q') in the combined n and q space. Once 
more, only significant eigenvectors of the norm matrix 
are retained. Finally, the Hamiltonian kernel is diago- 
nalized in the basis of norm eigenstates to construct the 
weight functions g). The final eigenstates \ JMi'), 

Eq. mix all K and q values. They are then used to 
calculate all observables and transition moments. 

Note that our method is not restricted to the choice 
of the triaxial quadrupole moment as a generator coor- 
dinate. This choice is suggested by the importance of 
quadrupole correlations in nuclei, but any other collec- 
tive variable associated with an even multipole moment 
could be additionally considered, at the expense of hav- 
ing a larger mean-field basis. Inncluding odd multipoles, 
in particular octupole deformations, would require a rel- 
atively simple generalization of the codes but cannot be 
considered at present. However, some modes not related 
to a shape degree of freedom, such as pairing vibrations, 
can be included without modification of the numerics, 
as done earlier without angular momentum projection in 
Refs. [ii,[2i,[33]. 

D. Electromagnetic matrix elements 

Once the HWG equation is solved and the weight func- 
tions F^{n,q) are known, the expectation values and 



transition moments of other observables O can be de- 
termined starting from the kernels {JMKq\0\JMK'q') 
of the corresponding operators. Some of them provide a 
good test of the accuracy of symmetry restoration, which 
can be used to determine a sufficient number of points 
for projection. The matrix elements of the proton and 
neutron number operators, for example, are equal to the 
required values with an absolute deviation lower than 
10~^, and they have a dispersion lower than 10~* when 
using a nine-point formula for the particle-number pro- 
jection. 

The symmetries of the unrotated time-reversal invari- 
ant HFB states (time reversal and triaxiality) permit the 
reduction of the integration intervals for Euler angles to 
a e [0,7r/2], /3 e [0,7r/2], and 7 e [0, tt] (i.e., 1/16 of 
the full Svr^ integration volume for systems with integer 
spin). The integrals over a and 7 are discretized with a 
trapezoidal rule, while for (3 wc employ Gauss-Legendre 
quadrature. The number of points in these intervals used 
for the results reported in the following is 6 for a, 18 for 
/3, and 12 for 7, which corresponds to 24, 36, 24 points 
in the full Stt^ integration volume necessary for integer J 
values. With this discretization, the calculated expecta- 
tion values of the angular momentum operators and 
Jz are accurate with an error of the order 10^"' for the 
values of J discussed here. 

Since the restoration of rotational symmetry provides 
wave functions in the laboratory frame of reference, one 
can compute directly the expectation values and transi- 
tion matrix elements of electromagnetic operators. Be- 
sides mean-square radii and EO transition moments, we 
calculate spectroscopic quadrupole moments and B{E2) 
values as well as magnetic moments and B{Ml) values. 
By combining projection and variational GCM mixing 
of states with different mass quadrupole moments, the 
B{E2) transition moments for in-band and out-of-band 
transitions take the form 

B{E2- 4, ^ J,) 
= TTTTT E E E \{JM.\Q,,\J'MV)\' 

M=-J M' = -J' n=-2 

(11) 

and the spectroscopic quadrupole moments are given by 



Qs{Ju) = Y — = J v\Q2q\JM = Jv). (12) 

The matrix elements of the electric quadrupole moment 
operator Q2fi = e X)p '^p ^2^(^^p) that enter the B{E2) 
values and Qs are calculated by using point protons with 
their bare electric charge e. This approach is justified by 
the following considerations. First, there is no empirical 
evidence that electric moments are significantly modified 
by the in-medium effects that are absorbed into the ef- 
fective energy functional but ignored for all other observ- 
ables. Second, we use the entire space of occupied single- 
particle states without assuming an inert core. Third, 
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the projected GCM, as a "horizontal expansion" [53| of 
large-ampUtude dynamics, can be formulated such that 
the calculation of matrix elements does not contain a sum 
over unoccupied states [13, [13, [EH] • The absence of effec- 
tive charges is an important feature of a method aiming 
at a universal description of nuclei. The (at least ap- 
proximate) folding with the intrinsic charge distribution 
of protons and neutrons to construct a charge density 
is important for radii but it plays no role for mul- 
tipole moments when accepting the current precision of 
our approach. 

Both the B{E2) and Qs values scale with mass and 
angular momentum. With triaxial shapes, however, on 
can no longer define a dimensionless intrinsic deformation 
P2 that corresponds to a single Qs or a single B{E2) 
value, as can be done in the axial case [35| . as one needs 
two independent observables (i.e., that are not within 
the same band) to determine the deformation (32 and the 
angle 7 [56|. 

Similarly, magnetic moments are given by 

^(J,) = (J, M = J, v\fiio\J, M J, u) (13) 
and B{M1) values are obtained as 

= E E E \{JM.\^.\J'M'.t 

M=-J M' = -J' 

(14) 

The magnetic dipole operator fj, entering both expres- 
sions is given by 



(15) 



where Lj and St are the total orbital and spin opera- 
tors for protons and neutrons, t — p,n, and gi^p = 1 ^n, 
gs,p = 5.585 fiN, and gs,n = — 3.826 /^at are the g factors 
of protons and neutrons in units of the nuclear magneton 
fiN- Unlike the electric moments, the spin contribution 
to the magnetic moments is modified by the physics of 
short-range correlations which is resummed into the en- 
ergy functional but not explicitly considered in the wave 
functions. The spin g factors are often quenched by an 
empirical factor to mock up these effects. In the present 
paper, we start with time-reversal invariant mean-field 
states for which the spins of conjugated states exactly 
cancel each other. Although this is no longer the case 
after mixing such states, particularly when projecting on 
angular momentum, the angular-momentum projected 
total spin {J M=J i/\Sz\J M—J i^), remains a small cor- 
rection of a few percent to the dominating collective con- 
tribution from orbital angular momentum when calculat- 
ing the z component of total angular momentum in the 
laboratory frame, {J M=J vl^lJ M=J = HJ. This 
has two consequences. First, we can ignore the subtlety 
of in-medium corrections to the spin g factors for time be- 
ing. However, when projecting and mixing cranked states 



at high spin or n-quasiparticle states that might have sig- 
nificant open spin, this point will deserve more attention. 
Second, the magnetic moments calculated for the pro- 
jected and mixed states discussed in the present paper 
will be very close to the expectation value of the proton 
orbital angular momentum {J M—J vlLz.plJ M—J i^) in 
nuclear magnetons. 



E. Evaluation of kernels 

The matrix elements that are needed to solve the 
HWG equation and to calculate the electromagnetic 
properties of the GCM eigenstates are matrix ele- 
ments of tensor operators Tx/i of rank A between 
particle-number and angular-momentum projected states 
{q\pNpzpj^f^^pr^^pzpN^^,y As at present we 

consider only operators that commute with the particle- 
number operator, the particle-number projector for the 
left state can be commuted with Txf^ and one only has to 
project the right state on particle number. The angular- 
momentum projection operator, however, does not com- 
mute with tensor operators of rank different from zero. 
Still, the commutator can be evaluated using angular- 
momentum algebra to obtain matrix elements with an 
angular-momentum projection operator acting on the left 
state only: 



2J' + 1 
27" 



+J 



- {J'XJ\M'nM) J2 {J'XJ\K',K' ~k,k) 



k=-J 



■{q\PLTx^K'-kW) 



(16) 



For a scalar operator Tqq, this simplifies to 

{qlPLlTooPlpK'W) = Sjj.SMM'iqlP^K' TooW) ■ 

(17) 

As a consequence, we have to evaluate matrix elements 
of the form 

{q\R\a, f3, 7) Tam e'*"^ e^*'^\q') = {L\fx^\R) (18) 

where \q) and \q') are HFB states with possibly differ- 
ent intrinsic deformations, and |L) and |R) are the HFB 
states rotated in coordinate and gauge space, respec- 
tively. Key elements of the evaluation of these kernels 
are outlined in the Appendix IIVI 



F. Configuration mixing using energy functionals 

For the sake of simplicity, we have outlined the pro- 
jected GCM using a many-body Hamilton operator. In 
practice, however, an energy density functional is used 
in a for m p ostulated by analogy with the Hamiltonian 
case [13,[53]- In this procedure, the multi-reference (MR) 
energy density functional (EDF) is obtained by replac- 
ing the density matrices in the single-reference energy 
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FIG. 1: (Color online) Lower panel: Nilsson diagram of the 
eigenvalues of the proton single-particle Hamiltonian for axial 
shapes. Solid lines denote levels with positive parity, dotted 
lines levels with negative parity, and the red dashed line the 
Fermi energy. Upper panel: Particle-number projected (solid 
line) and particle-number and angular-momentum J — pro- 
jected (dotted line) axial deformation energy curves. 



for the particle-hole channel of the effective interaction, 
supplemented by a zero-range density-dependent pairing 
energy functional [6^. A soft cutoff around the Fermi 
energy is used when solving the HFB equations as intro- 
duced in Ref. [G^I • 

Figure [T] provides as a reference some key results from 
an axial calculation as a function of the dimcnsionless 
mass quadrupole deformation given through Eq. ([3]). In 
this figure, oblate shapes are labeled by negative values 
of P2- On the top part we compare the variation of the 
energy with deformation obtained by projecting on parti- 
cle numbers N and Z only (solid line) and by projecting 
on particle number and on J = (dotted line). In the 
bottom part, the Nilsson diagram of eigenvalues of the 
proton single-particle Hamiltonian is shown . The Nils- 
son diagram for neutrons is nearly identical to the one for 
protons, except for a global shift owing to the absence of 
the Coulomb potential. The projection on J = signifi- 
cantly increases the energy gain from deformation. The 
mean-field configuration corresponding to the minimum 
also has a significantly larger prolate deformation after 
projection than before. After projection, the axial mini- 
mum corresponds to the intrinsic deformation where the 
Z — N = 12 gap in the Nilsson diagram is largest. The 
Nilsson diagram indicates also that there are only very 
few level crossings in the interesting region of deforma- 
tions. 



functional by the transition density matrices as they ap- 
pear when transition matrix elements of a Hamilton op- 
erator are evaluated with the generalized Wick theorem 
d, HO, Hi- The terms in the Skyrme EDF that con- 
tain time-odd densities are treated as described in Ap- 
pendix C of Ref. [1^ and give a small contribution to 
the nondiagonal kernels. The density-dependent terms 
are generalized by using the standard prescription that 
the density is replaced by the transition density in the 
density-dependent terms. This is the only prescription 
that guarantees various consistency requirements of the 
energy functional [H, [s^ . However, this procedure may 
lead in some cases to problems that have become evident 
recently (60l . [sij . Discontinuities or even divergences of 
the energy kernel may indeed appear as a function of 
deformation. A regularization scheme has been proposed 
62 1 and applied to a simple case of MR-EDF calculations 
63l |. However, it cannot be applied to the standard form 
of the Skyrme and Gogny interactions (g^l- We leave 
this problem to be addressed in the future and stick to a 
standard Skyrme energy functional, as done in our earlier 
studies. 



III. 2*MG AS AN ILLUSTRATION 

The main features of our method are exemplified by a 
calculation for ^^Mg. We use the parametrization SLy4 
[65| of the Skyrme-|-Coulomb energy density functional 



A. Triaxial energy maps 

Three {(3, 7) deformation energy surfaces for ^"^Mg are 
displayed in Fig. [21 All energies are normalized to the 
energy of the spherical configuration. A projection on 
particle numbers N and Z is performed for all three pan- 
els; it is combined with a projection on the angular mo- 
mentum J = in the middle panel and on J = 2 and 
K — on the lower one. Several comments on the inter- 
pretation of the deformation energy surfaces are in order. 
First, the coordinates (32 and 7 are not related to any ma- 
trix element in the angular-momentum projected energy 
surfaces but merely provide a label for the intrinsic state 
fromwhich the projected energy is obtained. All states 
projected on J = are spherical in the laboratory frame. 
Second, neither the mean-field nor the projected states 
are orthogonal, such that the energy maps do not and 
cannot represent the actual metric, which is related to 
the inverse of the respective overlap matrix. Finally, the 
J = 2, K = surface depends on the orientation chosen 
for the principal axes of the nucleus. The quadrupole 
moment along the z axis is the largest one in Fig. [51 We 
shall show later that the K = results have no obvious 
interpretation for other choices. 

Starting with the energy surface without angular- 
momentum projection, one sees that the mean-field 
ground state corresponds to an axial prolate deformation, 
more bound than the spherical configuration by about 
700 keV. Comparing with the axial energy curve given 
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FIG. 2: Triaxial quadrupole energy maps obtained by pro- 
jecting mean-field configurations on A'^ and Z only (top) and 
also either on J = (middle) or on J = 2 and A" = 0. In 
the latter case, the z axis is chosen as symmetry axis. Filled 
circles denote minima; the ® symbols denote saddle points. 
All energies are normalized with respect to the energy of the 
particle-number projected spherical state. 

as a solid line in the upper panel of Fig. [U one notices 
that the prolate minimum is indeed a true one but that 
the oblate extremum seen in Fig. [T] is a saddle point in 



the 7 direction. A qualitatively similar deformation en- 
ergy surface was found in HF-I-BCS calculations with the 
Skyrme interaction SIII in Ref. [sHi and in HFB calcula- 
tions with the Gogny force Dl in Ref. [695 . 

The projection on J = favors triaxial configurations: 
The lowest energy is obtained for a triaxial mean-field 
configuration with 7 w 16° and a value for fii around 
0.6, similar to that of the axial mean-field state giving 
the lowest J = projected energy. Before projection, 
these two states were separated by 2.2 MeV. The energy 
gain for the triaxial point is 3.1 MeV larger than that 
of the axial point. This difference is thus large enough 
to compensate for the difference in energy between the 
mean-field configurations. However, this result has to be 
taken cum grano salts as the deformation has a limited 
meaning after projection. In particular, the states result- 
ing from the projection on the same J value of mean- field 
configurations corresponding to different quadrupole mo- 
ments may have a large overlap. The states with the low- 
est energy obtained after projection of axial and triaxial 
configurations have an overlap close to 0.9. 

^''Mg is one of the few light nuclei with a deformed 
mean-field ground state when pairing correlations are 
taken into account. In this small system, however, the 
static quadrupole correlation energy (i.e., the deforma- 
tion energy of the mean-field ground state) is much 
smaller than the dynamical quadrupole correlation en- 
ergy obtained from projection and mixing of states with 
different intrinsic shapes; this is a general feature of light 
nuclei with v4 < 100 "70]. 

The projection on J = 2, X = leads to an energy 
map whose interpretation is more difficult than that for 
J = 0. The orientation of the nucleus that is chosen gives 
to this energy map some meaning in the interpretation 
of the first 2+ level, whereas the second 2+ state belongs 
mainly to the map corresponding io K = 2. The topog- 
raphy of the J = 2, K = Q map is quite similar to, but 
slightly flatter than, that of the J = energy map. 

We will not show here the deformation energy sur- 
faces for higher J and K values. When starting with 
time-reversal invariant mean-field states, the angular- 
momentum projected energy surfaces for all even values 
of J and K ^ Q are infinite at 7 = 0°, and degener- 
ate with the K — Q surface of the same J for 7 = 60°, 
which we will exemplify for one oblate configuration in 
the next section. As a consequence, K mixing becomes 
usually much stronger when going from prolate to oblate 
shapes with increasing 7. For odd values of J, there is 
no if = component, and the energy surfaces for finite 
K are infinite at both 7 = 0° and 7 = 60°. 

These findings are consistent with what has been 
seen before. Many even-even nuclei have a quadrupole- 
deformed mean-field ground state with a substantial en- 
ergy gain that might be as large as 25 MeV in very 
heavy nuclei [TOj . The corresponding shape in the in- 
trinsic frame, however, remains axial in most cases. Tri- 
axial ground-state deformation is rare at the mean-field 
level, and the additional energy gain from nonaxial de- 
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FIG. 3: Excitation spectra and selected B{E2) values after decomposition of the mean-field state that gives the minimum of 
the J = projected energy surface into angular-momentum projected components \ JMKq) (left) and after K mixing into 
\JMhiq) (right) for three different orientations of the same triaxial state in the intrinsic frame. 



formation is usually quite small compared to the defor- 
mation energy of the axial minimum or saddle in the 
same nucleus (gI, [7l|, Izl, III [zi| • Triaxiality, however, 
plays an important role for high-spin states when de- 
scribed in mean-field approaches [tsI . [tH] • The situation 
changes when going " beyond the mean field" . It is well 
known that projection on good quantum numbers favors 
the breaking of the corresponding symmetries in the un- 
derlying mean- field states: For a nucleus with a spheri- 
cal minimum in the mean-field deformation energy sur- 
face, angular-momentum projection on J = results in 
a deeper minimum at slightly deformed shapes i7Qj, as 
first noticed by Dalafi [771; ^ nucleus with an axial 
minimum of the mean-field energy surface, the minimum 
of the J = angular-momentum projected energy sur- 
face is slightly shifted into the triaxial plane, as has been 
demonstrated by Hayashi et al. [lll |. 



B. Angular- momentum projection of a single 
triaxial HFB state 

1. The role of the orientation in the intrinsic frame 

Scalar quantities such as the overlap and the energy 
kernels of the GCM do not depend on the orientation 
of the projected state in the laboratory frame (i.e., on 
its angular momentum projection M). Moreover, final 
results do not depend on the orientation of the mean-field 
state in its intrinsic frame when time-reversal invariance 
is preserved. This equivalence, however, is only obtained 
after K mixing. This provides the opportunity to analyze 
the role of K mixing for the projection of states with 
same intrinsic deformation but a different orientation of 
their principal axes. It also constitutes an excellent test 
of the numerical accuracy of the projection scheme that 
we have developed. 



There are six possible ways to orient a triaxial nucleus 
in its intrinsic frame in such a way that the major axes of 
the intrinsic quadrupole moment coincide with the axes 
of the intrinsic coordinate system. These six possibilities 
correspond to the six sextants of the /3-7 plane. Ow- 
ing to the particular role played by the intrinsic z axis 
in the angular-momentum projector results are in- 
variant under exchange of the x and y axes such that 



TABLE I: Excitation energies after decomposition of the 
mean-field state that gives the minimum of the J = pro- 
jected energy surface into angular-momentum projected com- 
ponents \JMKq) and after K mixing into \ JMnq) for three 
difi'erent orientations of the same triaxial state in the intrinsic 
frame, corresponding to Fig. [3] 
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FIG. 4: Decomposition of the mean-field state that gives the 
minimum of the J = projected energy surface into compo- 
nents of same angular momentum J, Eq. (I19|) . The results 
obtained from the different possible orientations in the intrin- 
sic frame are undistinguishable. 

they give pairwise the same decomposition \JMKq). It 
is, therefore, sufficient to consider the intrinsic z axis co- 
inciding with the longest, intermediate, or shortest axis 
of the triaxial state, without specifying the orientation 
of the other two axes. The results for excitation energies 
and B{E2) values obtained for the three possible orienta- 
tions of the nucleus are presented in Fig. [3l key numbers 
for excitation energies are repeated in Table HI When the 
intrinsic long axis is chosen along the z axis (left panel), 
K mixing has a small effect and a clear connection can be 
made between the K — and the K = 2 bands and the 
ground-state and 7 bands after mixing. The differences 
between the results before and after K mixing are due 
to the nonorthogonality of the states | JMKq) with the 
same J and M values but different K , which is removed 
in the orthogonal basis {JMnq). This orthogonalization 
pushes up the second 4+ state by 1 MeV. The dominat- 
ing in-band B{E2) values are similar before and after 
K mixing, whereas the much smaller out-of-band B{E2) 
values change substantially. 

The situation is quite different when the triaxial mean- 
field state is orientated in such a way that the z axis is 
not the longest one. The energies and transition proba- 
bilities obtained before K mixing have no obvious corre- 
spondence with the K-mixed results. In particular, the 
B{E2) values without K mixing are of similar size for in- 
band and out-of-band transitions. However, both Fig. 
and Table |T] clearly indicate the independence of the re- 
sults from the orientation of the mean-field state after K 
mixing. Note that the 3"*" level is not affected by K mix- 
ing: as with good time-reversal invariance, there is only 
one independent nonzero K component with K = 2. 

Figure [3] presents the decomposition of the norm of 
the same states as in Fig. [3] in components with different 
total angular momenta J, 

4= E ^i\PiK\<i) (19) 

K=-J 

summed over all possible values of K\ Fig. [5] shows the 



FIG. 5: Decomposition of the mean-field state that gives the 
minimum of the J = Q projected energy surface into compo- 
nents of same K, Eq. (|20|) . for three different orientations of 
the state in the intrinsic frame. 

decomposition into components with different intrinsic 
angular momentum projection 

4= E ^^\Pkk\<i) (20) 

J>\K\ 

summed over all possible values of J for a given value of 
K. The underlying state |g) is assumed to be particle- 
number projected and normalized. As we always choose 
an orientation of the triaxial state \q) where the intrinsic 
z axis coincides with one of the principal axes, the com- 
ponents of opposite K have the same weight c^^^ = c^^. 
Furthermore, all components with an odd value of K are 
zero for all values of J. 

The decomposition of a triaxial mean-field state into 
components with different J values in the laboratory 
frame should be independent of its orientation in the in- 
trinsic frame, which is indeed the case within the reso- 
lution of Fig. m The plot suggests a separation of the 
coefficient c^j into two distinct curves: one for even val- 
ues of J, which dominates the decomposition and peaks 
for J — 2, and a second much weaker one for odd values 
of J, which peaks at J = 3. (When decomposing a time- 
reversal invariant state, there is no component with J = 1 
for symmetry reasons.) Increasing the deformation of the 
intrinsic state modifies the Cj distribution such that the 
peaks of the distributions for even and odd J are shifted 
toward larger values of J. 

By contrast, the decomposition of the same triaxial in- 
trinsic state into its K components depends strongly on 
its orientation, which underlines that K is not an ob- 
servable quantity. The distribution of K components is 
quite narrow when the long axis of the intrinsic state co- 
incides with the z axis, and becomes broader when the 
intermediate or even short axis is chosen instead. The 
different distributions in K space indicate that the nu- 
merical convergence of angular-momentum projection is 
not the same for all possible orientations of the mean- 
field state. The accurate determination of components 
with large K values requires many integration points for 
the Euler angles a and 7; hence, it is more favorable to 
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FIG. 6: Excitation spectra and selected B{E2) values af- 
ter decomposition of the mean-field state that gives the ax- 
ial oblate saddle of the J = Q projected energy surface into 
angular- momentum projected components \ JMKq). The left 
panel provides the decomposition when choosing the z axis 
as the symmetry axis, 7 = 180°. The middle panel gives the 
decomposition into K components when the symmetry axis 
is chosen perpendicular to the z axis. The right panel shows 
the unique band resulting from K mixing (see text). 



TABLE II: Spectroscopic quadrupole moments obtained from 
projection without and with K mixing compared with values 
obtained from the asymmetric rotor model. In this latter case, 
the intrinsic charge quadrupole moment Qo,p = 63.45 e fm^ 
(/3p = 0.583) and 7 = 16.1° value of the triaxial HFB state 
are used as input. 



Quantity 



No mixing K mixing Rotor 



Qs (2+) (e fm^) 




-19.1 


-18.7 


-17.0 


Qs (2+) (e fm^) 




20.0 


18.1 


17.0 


B(S2,2+ ->0+) 


(e^ fm^) 


94.9 


95.7 


75.5 


B(S2,2+ ^ 0+) 


(e^ fm^) 


15.9 


5.2 


4.6 


B{E2,2+ ^2+) 


(e^ fm*) 


8.4 


26.4 


14.3 


B{E2,3+ ^2+) 


(e^ fm*) 


11.6 


10.0 


8.3 


B{E2,3t ^2+) 


(e^ fm*) 


146.4 


165.4 


134.7 



one ends up with the same rotational band as obtained 
by exploiting the symmetry of the intrinsic state. 



3. Comparison with the asymmetric rotor model 



orient the mean-field wave function with its z axis along 
the long axis. The sum over all Cj coefficients should 
be equal to the sum over all the cj^ ones and equal to 
1. There is, in practice, a slight numerical deviation of 
the order of 10~^, which can be attributed to high-J and 
high-i^ components requiring a higher number of inte- 
gration points than have been used here. 



2. Decomposition of an oblate HFB state 

A further test of our method is given by the projection 
of an oblate mean-field configuration, for which the sym- 
metries of our codes allow two distinct orientations in its 
intrinsic frame: The z axis can be cither the symmetry 
axis or perpendicular to it. Using the properties of the 
transformation of the operator P//^ under rotation, one 
can show that a pure K = state is transformed into a 
multiplet of states with K between and J. However, 
the transformed wave functions differ only by an unob- 
servablc phase and are degenerate. An example is given 
in Fig. [5] for the mean-field state with mass quadrupole 
moment Qo — 80 fm^ which gives the axial oblate saddle 
point of the J = projected energy surface. When using 
the z axis as symmetry axis, as was done in axial calcula- 
tions tM,t35i,^, Ji?., 38i,|3l,|40,|4l[, the projection decom- 
poses the mean-field state into a single rotational band 
of K = states; all other K components have zero norm. 
When orienting the intrinsic state such that the z axis is 
perpendicular to the symmetry axis, as done at 7 = 60° 
in the (3-j plane, all even K components \JMKq) up to 
K ~ J are nonzero for each J . When constructing the 
\JMKq) states, however, the diagonalization of the norm 
matrix gives only one nonzero eigenvalue per even J, and 



The spectroscopic quadrupole moments and B{E2) 
values obtained by the projection of a triaxial mean- 
field state are mainly determined by its geometry. This 
property is illustrated in Table |TT] where we compare the 
values obtained by projecting the same triaxial mean- 
field state as in Fig. [3] to those calculated with the 
asymmetric rotor model introduced by Davydov and co- 
workers [tsI . [7^, using the intrinsic proton quadrupole 
moment Qo,p = 63.45 e fm^ (/3p = 0.583) and 7 = 16.1° 
as input. The agreement is excellent and, in practice, 
improves with deformation. 

C. Configuration-mixing calculations 

1. Selection of the mean-field basis 

As a last step, we perform a mixing of states ob- 
tained by projecting on particle-number and angular- 
momentum mean-field wave functions covering the full 
/3-7 plane. Specifically, the results will be analyzed by 
comparing the spectra and transition probabilities ob- 
tained in calculations using four different subspaces of 
states 

• a basis of "prolate" axial states, comprising the four 
deformations (91,92) = (80,0), (120,0), (160,0) and 
(200,0) fm2; 

• a basis of "axial" states, where two oblate deforma- 
tions (91,92) = (0,80) and (0,120) fm^ are added 
to the prolate basis; 

• a basis of six "triaxial" mean-field configurations 
(91,92) = (80,40), (120,40), (160,40), (200,40), 
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(80,80) and (120,80) fm^; 

• a basis labeled "full" where, depending on J, two 
to four states of the prolate basis are added to the 
triaxial basis. 

For the full basis, we have added for each J value the 
largest possible number of axial points to the triaxial ba- 
sis. We have excluded from the GCM calculations those 
deformed states that are situated in a region affected by 
the difficulties mentioned in Sec. Ill Fl This region spreads 
from small deformations around the spherical point to a 
region with larger deformations between the oblate axis 
and 7 « 50°. This restriction does not permit the mixing 
of triaxial and oblate axial states. 

For ^'*Mg, the number of axial states that can be added 
to the set of triaxial states is not large, just two to four, 
depending on the value for J. Redundancies appear very 
quickly in the norm kernel when more states are added 
to the nonorthogonal basis. This very small number of 
nonredundant states is a direct consequence of the very 
few level crossings visible in Fig.[Tl even states with quite 
different intrinsic deformation might have a very similar 
single-particle spectrum. This feature is illustrated in 
Fig. [71 where the eigenvalues of the J — norm kernel 
are plotted as a function of the number of states included 
in the configuration-mixing calculation. We start on the 
left of the figure with the eigenvalues of calculations per- 
formed in the prolate and axial bases. The addition of 
oblate configurations to the prolate set brings one large 
eigenvalue, close to 1, and another one around 10~^. The 
range of values obtained for a purely triaxial set is very 
similar to the axial set, although both bases have no vec- 
tors in common. Starting from the initial set of triax- 
ial points, prolate points are added one after the other. 
There is, thus, one more eigenvalue after each addition 
and the trace of the norm, which is equal to the number of 
discretization points, is increased by 1. One can see that 
this increase of the trace is mainly distributed among the 
largest eigenvalues, which increase slightly. Each time, a 
new eigenvalue around 5 x 10~^ also appears, or even an 
smaller one when the fourth axial point is added. Al- 





FIG. 7: Variation of the eigenvalues of the J — norm kernel 
among different bases. 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 

FIG. 8: (Color online) Deformation energy curve projected on 
and Z (black) only and projected also on J = (red) for 
three different values of 7, as indicated. The energies of the 
J = GCM ground state obtained from the axial, triaxial, 
and full bases as described in the text cannot be distinguished 
within the resolution of the plot and are represented by the 
same filled diamond plotted at arbitrary deformation. AH en- 
ergies are normalized to that of the spherical particle-number 
projected state. 



though these changes are very small, the effect of the 
coupling between axial and triaxial states on all observ- 
ables is not completely negligible. 



2. Ground-state correlation energy 

We have shown in Fig. [2] that the angular-momentum 
projection changes the topography of the J = defor- 
mation energy surface and generates a minimum for the 
projection of a triaxial configuration. This result is il- 
lustrated further in Figure [5] where the variation of the 
energy along three cuts in the plane is plotted. The 
first curve corresponds to prolate deformations, 7 = 0°, 
the second to oblate ones, 7 — 60°, and the third to a cut 
along 7 = 16°. Results obtained with and without pro- 
jection on J = are given. The big diamond covers the 
range of the GCM energies obtained by using the axial, 
triaxial, and full bases. The lowest energy corresponds 
to a prolate configuration in the nonprojected calcula- 
tion, with triaxial energies always much larger than the 
prolate and oblate ones. After projection on J = 0, the 
triaxial curve is below the prolate one for a large range of 
deformations and the absolute minimum corresponds to 
a triaxial configuration about 800 keV more bound than 
the axial saddle point. However, the total energy gain 
obtained by mixing axial configurations is larger than 
that from the projection of a single triaxial configuration. 
Moreover, there is only a 35 keV difference between the 
energies obtained by the mixing of triaxial configurations 
and a further gain of 160 keV in the largest possible set 
of axial and triaxial configurations. 

This result puts a limit on the meaning of " rotational" 
and "vibrational" energies in the ground state: What 
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FIG. 9: Excitation spectra and B{E2) values in fm'* obtained in three of the bases defined in the text compared to the 
available experimental data. In the four cases, the spectrum is subdivided into a ground-state band, a 7 band, and additional 
low-lying states that do not necessarily form a band. 



appears to be "vibrational" energy in a projected ax- 
ial quadrupole GCM is transformed into " rotational" en- 
ergy in a projected GCM including triaxial states. At 
the same time, this result is rather encouraging as it 
brings the hope that the ground-state correlation ener- 
gies are rapidly saturating when more than the restora- 
tion of symmetries on the axial mode is included. Hagino 
et al. [soi l arrived at a similar conclusion while studying 
the ground-state correlation energy in an exactly solv- 
able model that simulates collective quadrupole motion. 
Of course, systematic realistic calculations, in particular 
for heavier nuclei, are needed before a general rule can 
be established. 



3. Excitation spectrum 

The excitation spectra obtained in three different 
GCM bases are given in Fig. M and Table [ml Table |TTT] 
also lists results for spectroscopic quadrupole moments, 
while the B{E2) values for the most important transi- 
tions are given in Fig. [3) Experimental data are taken 
from Ref. f8ll| in both cases, using 5.940 x 10"^ ^4/3 g 
fm"' = 4.11 e fm''^ for the E2 Weisskopf unit. 

The spacing of levels, the relative strength of B{E2) 
values and the K content of the collective wave functions 
suggest a separation of the low-lying spectrum into three 
bands: The AJ = 2 ground-state band dominated by 
K = components, a A J = 17 band dominated by K — 
2 components, and a third band that is again dominated 
by if = components. 

As already mentioned, the GCM ground-state energy 
is quite close in the purely axial and purely triaxial cal- 
culations, but is 160 keV lower in the full calculation that 
combines triaxial and axial states. All excited levels ex- 
cept the O2 one, are lower in the triaxial basis than in 
the axial one. Adding prolate axial states to the triaxial 



basis mainly lowers the energy of all 0+ levels, push- 
ing up levels with other values of J. The spectroscopic 
quadrupole moments Qs and B{E2) values, however, are 
not significantly affected. 

The most significant difference between the axial and 
triaxial calculations concerns the first excited band, 
which is clearly a 7 band in the triaxial basis. Odd-J 
levels are of course absent from the first excited band in 
the axial basis but some of the features of a 7 band are 
already hinted in this band dominated by oblate configu- 
rations. In particular, the B{E2) value for the transition 
from the 4^ to the 2^ level is very close to that of the 
triaxial calculation. This identification of the projection 
of an oblate state to a 7 band was already suggested 
in Sec. IIII B 21 However, K mixing does not bring any 
gain in energy in the projection of an oblate state since 
the relative weights of the K components are fixed by 
symmetry relations, whereas they are free to vary in the 
projection of triaxial states. The gain in energy brought 
by the inclusion of triaxial configurations is evident when 



TABLE III: Comparison between the theoretical excitation 
energies in MeV and spectroscopic quadrupole moments in 
units of e fm^ and the experimental values taken from 
Ref. B. 

Level Eex Qs 





Axial Triaxial 


Full 


Expt. 


Axial 


Triaxil 


Full 


Expt. 




2.24 


1.87 


1.97 


1.37 


-17.1 


-19.6 


-19.4 


-16.6 (6) 




5.88 


5.44 


5.57 


4.12 


-25.1 


-26.1 


-26.0 






7.69 


6.88 


6.99 


4.24 


9.9 


17.1 


16.6 




it 




9.59 


9.74 


5.24 




-0.1 


-0.1 






13.29 


11.12 


11.28 


6.01 


9.0 


-7.3 


-7.4 






7.53 


8.79 


7.520 


6.42 


0.0 


0.0 


0.0 


0.0 
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TABLE IV: Overlaps between the GCM wave functions ob- 
tained in calculations using different sets of mean-field states. 
The two results marked by asterisks indicate that the number 
given represents the overlap between the second prolate and 
the third triaxial or full 2^ excited states. 







Triaxial 


Full 






Prolate Axial 


Prolate Axial 






0.97 


0.99 


0.97 


0.99 


J 


= 


0.89 


0.94 


0.96 


0.97 






0.16 


0.93 




0.96 






0.97 


0.98 


0.97 


0.98 


J 


= 2 


0.21 


0.86 


0.23 


0.88 






0.88* 


0.82 


0.88* 


0.81 



looking at the energies and the spectroscopic quadrupole 
moment f Table Hill) , which change drastically when going 
from the axial to the triaxial calculation. 

The inclusion of triaxial deformations has also a strong 
effect on the energies of the levels in the second band. 
The excitation energy of the 2 J level is lowered by nearly 
1 MeV, and that of the excitation energy of the 4^ level 
even more. Altogether, the inclusion of triaxial mean- 
field states brings the right tendency to make the spec- 
trum more compact and brings it closer to the exper- 
imental data. It remains to verify whether the use of 
time-reversal invariance breaking HFB states will add 
the extra gain of energy required to bring the states with 
J different from closer to the experimental levels. In 
particular, the odd-J states in the 7 band still have ex- 
citation energies that are much too high. 

To analyze further the equivalence of and differences 
between the bases that we have used, the component of 
the triaxial and full bases that are included in the pro- 
late and axial bases are given in Table [TVl These over- 
laps, (JMi'lJMfj,), can be easily calculated by using the 
FiJ{K,q) and the norm kernels T{K,q; k' ,q'). A bit sur- 
prisingly, the differences between the collective ground- 
state wave functions obtained within the axial and full 
bases are very small, the overlap with the prolate basis 
being slightly lower but still quite high. The second 0+ 
state is not as well described by the prolate basis and 
the third one state is missing in this basis. Note that the 



TABLE V: Comparison between theoretical and experimental 
B{M1) values in units of filf. Data taken from Ref. [sij . 



transition 


B{M1) 




Axial 


Triaxial Full 


Expt. 


2+^2+ 4 X 10"^ 

3+^2+ - 

3+^2+ — 
■^1 ^ ^2 


3 X 10"^ 3 X 10"*^ 
2 X 10"^ 1 X 10"'^ 

4 X 10"^ 4 X 10"^ 


1.6 X 10"^ (14) 
3.8 X 10"^ (20) 
6.2 X 10"* (30) 



difference between the axial and full ground-state wave 
functions is still larger than the energy differences be- 
tween these states (0.16 MeV out of a total energy of 
around 200 MeV). 

The first 2+ state has the same structure in all bases. 
In contrast, the second 2+ state of the axial and full 
bases is not described by the prolate basis. Although the 
overlap between the 2^ states in the axial and full bases 
is close to 0.9, the excitation energy is lowered by 800 keV 
when triaxial configurations are included. This confirms 
our previous interpretation that the mixing of prolate and 
oblate axial states can, to some extent, simulate states 
with K 0, but not fully. 

In the N = Z nucleus ^'*Mg, all calculated magnetic 
moments are just a few percent larger than ZJ/{N + 
Z) = J/2 in nuclear magnetons. This is a consequence 
of (a) the time-reversal symmetry that we impose on the 
underlying HFB states as explained in Sec. lIIDl and of (b) 
the fact that protons and neutrons have nearly the same 
contribution to the angular momentum in this N = Z 
nucleus. The calculated magnetic moments agree well 
with the available experimental ones [sij for the 2]^, 4^, 
2^, and 4^ states within the experimental error bars. 

The calculated and experimental B{M1) values are 
compared in Table [V] Data are taken from Ref. (sij 
where 1.790 is used for the Ml Weisskopf unit. The 
values that we obtain are about one order of magnitude 
too small, which clearly indicates that the projected cur- 
rents and spin densities are underestimated. Starting 
from time-reversal breaking mean-field states instead of 
time-reversal invariant ones as done here can be expected 
to increase the B{M1) values. 

IV. SUMMARY AND OUTLOOK 

A generalization of a method that enables the mixing 
of projected mean-field states that was previously lim- 
ited to axial configurations has been set up to allow for a 
description of the full five-dimension quadrupole dynam- 
ics. Compared to a GCM calculation limited to projected 
axial quadrupole deformed states, the present method al- 
lows for a spreading of the collective states into the (3-^ 
plane. In the case of '^'^Mg, we have shown how the spec- 
troscopic properties of the low-lying states are affected. 
One can summarize our main results as follows: 

• When looking at the projection of a single config- 
uration, the energy obtained for the ground state 
is significantly lowered when allowing for triaxial 
quadrupole deformation; 

• If one considers at the same time the correlation 
energy from symmetry restorations and configura- 
tion mixing, the total energy difference between the 
ground state obtained within an axial and a triaxial 
mean-field basis is quite small. This indicates that 
the nondiagonal matrix elements between prolate 
and oblate axially deformed mean-field states bring 
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a large fraction of the correlation energy that is ob- 
tained by the projection of triaxial configurations 
between them. This puts a limit on the meaning 
of rotational and vibrational energies, as what ap- 
pears to be "vibrational" energy in the projected 
quadrupole GCM of axial states is transformed into 
"rotational" energy in a projected GCM including 
triaxial states. Of course, this result will have to 
be confirmed by studies in other nuclei; 

• This finding is also supported by the analysis of the 
spectrum of eigenvalues of the norm matrix, where 
adding triaxial states does not introduce states with 
large eigenvalues to the space of 0+ states. 

• For higher J values, the situation is more complex. 
However, including triaxial deformations lowers the 
excitation energies and brings the spectrum closer 
to experiment. 

the basis of this analysis, one can draw some conclu- 
about the effect of triaxial deformations: 

• Their effect on masses seems marginal and it is 
reassuring in some way: If one can confirm that 
triaxial deformations increase binding energies by 
only around 100 to 200 keV, it would be justified 
to avoid the complexities from their introduction in 
systematic mass calculations. 

• The gain of energy on excited states partly cures 
a problem common to all projected GCM calcula- 
tions based exclusively on axial mean-field states. 
However, further improvements are still necessary, 
including the breaking of time-reversal invariance 
and the consideration of the projection of crank- 
ing states optimized for each angular momentum. 
The breaking of axiality is a necessary first step 
before the breaking of time-reversal invariance and 
we hope to validate the extension of our method in 
this case in the near future. 

• The power of our method will be more apparent 
when breaking of time-reversal invariance is in- 
cluded. Although some new states can already be 
described at the present level of development (e.g., 
odd- J members oi K ^ Q bands), it will be possi- 
ble to describe quasiparticle excitations and in par- 
ticular nuclei with an odd number of neutrons or 
protons. 

• Our method also provides an ideal tool to bench- 
mark simpler models. To give only one example, 
the metric of the /3-7-plane is generated directly 
in our method by the overlaps between mean-field 
wave functions of different shapes. Although it is 
not trivial to derive this metric in a multidimen- 
sional problem, as was done in Ref. [l^l for the one- 
dimensional case, it would be very instructive to 
compare a metric derived from a purely microscopic 



approach to the metrics that are usually based on 
semiclassical approximations. 

There are no basic reasons that prevent the application of 
the method presented here to heavy nuclei, although this 
would be too time consuming with our current numeri- 
cal implementation. Further developments are required 
to improve the efficiency of the numerical algorithms. In 
particular, the choice of discretization points used to eval- 
uate the integrals over the Euler angles needs to be opti- 
mized and the codes have to be parallelized, which could 
be done in a very efficient way. Work in that direction is 
underway. 

Another development that will have to be completed in 
the near future is the implementation of the regulariza- 
tion scheme proposed in Ref. [6^ to remove the patholo- 
gies brought by the use of the generalized Wick theorem 
when evaluating the energy kernels starting from an en- 
ergy density functional. However, the conceptual and 
technical difficulties encountered in the present general- 
ization of the projected GCM justify the continuation of 
the present developments in parallel. Considerable work 
at the level of the computational algorithm still remains 
to be performed to have a method applicable to heavy 
nuclei. 
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APPENDIX: COMMENTS ON THE NUMERICAL 
EVALUATION OF THE KERNELS 

1. Self-consistent mean- field calculations 

The mean-field states are generated with a simplified 
version of the cranked HFB code whose evolution is doc- 
umented in Refs. [48j, [661, [68|. We have imposed time- 
reversal invariance of the HFB vacua such that the single- 
particle wave functions are pairwise connected by time re- 
versal. The single-particle wave functions are represented 
as complex spinors discretized on a three-dimensional 
Cartesian Lagrange mesh in coordinate space. In ad- 
dition to time reversal, two further symmetries of the 
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symmetry group [H, HI] are imposed on the single- 
particle basis, namely that they are eigenstates of the 
parity P and z signature Rz. Their relative phases are 
fixed by choosing a basis where the eigenvalue of the 
y time simplex Sy is -1-1 for all single-particle states. 
This introduces three plane symmetries and allows us to 
restrict the numerical representation of individual wave 
functions to 1/8 of a full box ^] . In addition, the single- 
particle states are chosen to have good isospin projection 
(i.e., they are either pure proton or neutron states). 

Mean-field states with different deformation are ob- 
tained adding constraints on qi and 52 as defined in 
Eq. ([T|) in the variation. The constrained HFB equations 
are solved by using the "two-basis method" described in 
Ref. [HI, which delivers the HFB states \q) through quite 
a small number of single-particle states represented in co- 
ordinate space in the so-called Hartree-Fock (HF) basis 
that diagonalizes the single-particle Hamiltonian and the 
corresponding U and V matrices that establish the gen- 
eral Bogoliubov transformation from the HF ba- 
sis to the quasiparticle basis that diagonalizes the quasi- 
particle Hamiltonian. This procedure permits to limit 
the numerical representation to all single-particle levels 
below the Fermi energy and to a small number of levels 
above. As done in our earlier configuration-mixing cal- 
culations, we add the Lipkin-Nogami (LN) prescription 
to the HFB equations to enforce pairing correlations in 
all mean-field states. Using states without pairing in a 
GCM calculation introduces the danger of artificially de- 
coupling many-body states with a different ratio of occu- 
pied single-particle states of positive and negative parity, 
which can lead to instabilities when solving the HWG 
equation ([8]). 

The representation of the single-particle states on a 
coordinate space mesh has the clear advantage that its 
precision is fairly independent of the deformation when 
a sufficiently large box is chosen, and it only depends on 
the distance of discretization points. This is important 
for GCM calculations mixing states with very different 
deformation. 

For the subsequent projection and mixing of HFB 
states with different shapes, it is advantageous to use the 
canonical single-particle bases of the mean-field states as 
starting point. This simplifies the corresponding U and 
V matrices and allows for a safe cutoff of single-particle 
states with very small occupation v'^ that do not con- 
tribute to any of the kernels. For the identification of 
symmetries of the integrals over Euler angles, it also turns 
out to be advantageous to transform the single-particle 
basis to a basis of eigenstates of x signature Rx- 



2. Rotation of mean-field states 

The rotation of the " left" state in coordinate space and 
the "right" state in gauge space, Eq. can be per- 

formed either as rotations of the canonical single-particle 
states leaving the corresponding U and V matrices un- 



touched or as rotations of the U and V matrices leav- 
ing the single-particle states untouched. For the coor- 
dinate space rotation, the latter is equivalent to the ex- 
pansion of the rotated single-particle states in terms of 
the unrotated ones, which already at moderate deforma- 
tion requires highly excited single-particle states above 
the Fermi energy, which are outside of the single-particle 
basis used to describe the unrotated state. Therefore, 
the coordinate space rotation R{a, P, 7) is performed as 
a rotation of the single-particle states on the mesh as 
described in Refs. [ij, [8J|, which, however, is the most 
time-consuming piece of the numerical calculations. In 
contrast, it is simpler and numerically more efficient to 
perform the gauge-space rotation as a transformation of 
the V matrices instead of the single-particle basis, which 
in the canonical basis boils down to the multiplication of 
a small antidiagonal matrix with a phase factor. 

The rotation operator R(a, /?, 7) mixes single-particle 
states of both signatures, which requires one to extend 
the numerical representation of the single-particle wave 
functions from 1/8 to 1/2 of the full box, leaving only 
parity (and the isospin projection) as a good quantum 
number. The symmetries of the unrotated time-reversal 
invariant HFB states, however, permit the reduction of 
the integration intervals for Euler angles to a S [0,7r/2], 
/3 e [0,7r/2], and 7 £ [0,7r] (i.e., 1/16 of the fuU Sn"^ 
integration volume for systems with integer spin). 



3. Calculation of the kernels between rotated 
mean-field states 

Rotating an HFB state in coordinate or gauge space 
gives back an HFB state; hence, the matrix elements be- 
tween "left" and "right" states, Eq. can be easily 
evaluated with the generalized Wick theorem [sl] . The 
expressions given in Refs. [I,!!!] cannot be directly used, 
however. First, we have to transform the contractions 
between states in different quasiparticle bases to expres- 
sions for contraction between states in different canoni- 
cal single-particle bases [20l. [22j. Second, our coordinate 
space representation has as a consequence that the single- 
particle bases that set up \q) and \q') are not equivalent, 
a difficulty that is amplified further by rotation of one of 
the states. This difficulty can be overcome by eliminat- 
ing the contribution to the kernels coming from single- 
particle states that are occupied in one of the bases but 
not in the other ^ 0, [2I, HI, HI] . 

Only diagonal kernels and half of the off-diagonal ones 
have to be calculated; the other half of the off-diagonal 



^ The expressions for the matrices defining the basic contractions 
and the overlap given in Refs. [20I . |23 . ISi contain a systematical 
typographical error: all ("Rt)-! should be replaced by (TZ'^)-^, 
which does not have any consequences for these papers as the 
matrix TZ is real with the symmetries chosen there. 
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kernels can be constructed by using symmetries of the 
kernels. 



4. Particle-number projection and the phase of the 
overlap 

The integrals over the gauge angles for projection on 
proton and neutron number are discretized with the 
Fomenko prescription [85| , which is equivalent to a trape- 
zoidal rule. By using the number parity of the mean- 
field states, the integration intervals can be reduced to 
4> G [0,7r] for protons and neutrons. Additionally, a sym- 
metry connects the basic contractions and overlap in the 
interval [0, 7r/2] with those in [7r/2,7r], which, however, 
cannot be used to reduce the integration intervals for the 
reasons outlined in the Appendix of Ref. [1^ . 



With the symmetries of the HFB states chosen here, 
the overlap of two HFB states with different deformation, 
where one is additionally rotated in coordinate space, re- 
mains real. As soon as one of the two HFB states is 
additionally rotated in gauge space, however, the overlap 
in general becomes coniplex. Its modulus is determined 
by the Onishi formula [s, To^l. Its phase, a rapidly vary- 
ing oscillatory function of the gauge angle of particle- 
number projection, has to be determined from continu- 
ity arguments. Its value is zero for the overlap between 
non-gauge-space- rotated states and is followed during the 
gauge-space rotation by performing a second-order Tay- 
lor expansion. To this end, the overlap and its derivatives 
are also determined at a small number of gauge angles 
between the integration points for the gauge-space inte- 
gration, separately for protons and neutrons. 
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